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The recently developed technique combining the weak-coupling limit with the Floquet formalism 
is applied to a model of a two-level atom driven by a strong laser field and weakly coupled to 
heat baths. First, the case of a single electromagnetic bath at zero temperature is discussed and 
the formula for resonance fluorescence is derived. The expression describes the well-known Mollow 
triplet, but its details differ from the standard ones based on additional simplifying assumptions. 
The second example describes the case of two thermal reservoirs: an electromagnetic one at finite 
C*") ' temperature and the second dephasing one, which can be realized as a phononic or buffer gas 

reservoir. It is shown using the developed thermodynamical approach that the latter system can 
work in two regimes depending on the detuning sign: a heat pump transporting heat from the 
dephasing reservoir to an electromagnetic bath or heating both, always at the expense of work 
f~\ , supplied by the laser field. 

<L> 

in _ 

The model of a two-level system (TLS) driven by a classical electromagnetic field - laser radiation and subject to 
dissipation/decoherence effects is one of the most studied models in quantum optics and related fields. The basic 
I . physical phenomena - Rabi oscillations and the related appearance of Mollow triplet in the fluorescence spectrum - 
can be explained by the simplest version of such models It is based on the Markovian master equation for the 
& i reduced density matrix of the TLS which could be written as (h = ks = 1) 

^± = -i[H(t),p(t)} + £p(t) (1) 
' where H (t) is the Hamiltonian of the TLS including the driving harmonic field 

m 
in 



I. INTRODUCTION 



H(t) = ^ <7 3 + g{a~e int + a + e~ int ) (2) 
and C the phenomenological generator of Lindblad-Gorini-Kossakowski-Sudarshan (LGKS) form 



Cp = ^([a-,pa+] + [a-p,a + ]) + ^([cr+ pa"] + [v+p,*-]) S -[a s , [a 3 ,p]} (3) 

with phenomenological damping (74,), pumping (7-f), and dephasing (5) rates. 

It long has been noticed that the above equation is not consistent with the second law of thermodynamics and also 
does not reproduce experimental results well for the strong driving case (seeRef. for the discussion and numerous 
references). The reason is that the dissipation and decohcrence mechanisms are not independent of the Hamiltonian 
dynamics, which is now time dependent and, hence, one should expect that the simple phenomenological generator ([3]) 
should be replaced by a time-dependent superoperator. There exist several examples in the literature addressing this 
problem for different physical systems like atoms, NMR, quantum dots, and superconducting qubits. The authors 
derive master equations using various approximation schemes 0]. In the present paper we use a new systematic 
approach which combines the Davies idea of the weak-coupling limit Q with the Floquet theory of periodically driven 
quantum systems. This approach has been outlined in Ref. Q and used for different models of quantum engines 
and refrigerators @ . As already noted, the mathematical structure of the obtained (time-inhomogeneous) Markovian 
master equation is consistent with the second-law of thermodynamics and allows us to define unambiguously heat 
currents and stationary power when applied to models of quantum machines. 

The motivation to improve the existing formalisms, as well as for the TLS dynamics governed by the Hamiltonian ([5]) 
(which has not been studied within this new framework so far), is as follows: 

(1) The experiments concerning control of microscopic quantum objects like single atoms, quantum dots, or nanoscopic 
superconducting devices become so precise that the commonly used approximations may be insufficient. 

(2) Recently, thermodynamics of microscopic quantum machines has been studied intensively and it is important that 
the theoretical formalism does not introduce spurious inconsistencies with thermodynamics which are entirely due to 
the chosen approximation schemes (see Ref. [7[ for a sample of references) . 
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(3) The advocated formalism produces rather simple canonical quantum master equations, which reduce to simple 
kinetic equations for populations in a proper orthonormal basis. This allows study analytically various new designs 
of microscopic quantum machines. 

We do not repeat here the derivation of the master equation for periodically driven quantum open systems, which 
can be found in Ref. @, but present only its main ingredients for its construction and basic properties without 
proofs (Sees. II and III). We then discuss the case of TLS with the Hamiltonian ([2]) applied to two particular 
models. The first one corresponds to the resonance fluorescence phenomenon, including detuning effects. Although 
the result qualitatively reproduces the well-known Mollow triplet, there are quantitative differences in comparison 
with the existing approximation schemes which, in principle, should be detectable. The second model presents a 
microscopic heat pump based on the driven TLS (atom, quantum dot, nanoparticle, etc.) which illustrates a new 
coherent mechanism of energy transfer among TLS, an electromagnetic field, and a generic dephasing bath, which 
involves Rabi oscillations. 



II. MASTER EQUATION FOR PERIODICALLY DRIVEN OPEN SYSTEM 

A "small" quantum system (i.e., described by a Hamiltonian with a discrete spectrum) is driven by an external 
periodic field and weakly coupled to an environment. The environment is assumed to be a "large" quantum system 
which can be treated as a system with continuous spectrum satisfying all standard ergodic conditions. They allow us 
to apply the weak-coupling limit procedure to derive the Markovian master equation for the reduced density matrix of 
the small system. The reversible effect of driving combined with the Hamiltonian corrections caused by the interaction 
with environment ("Lamb shift") is described by the physical (renormalized) periodic Hamiltonian H{t) = H(t + r). 
We assume the following form of the system-environment coupling: 

H int = S®F , (F) E = Q. (4) 

where the environment is assumed to be at the stationary state and (■ • - )e denotes the associated quantum average. In 
particular, in the weak-coupling Markovian approximation the influence of the environment is completely characterized 
by its spectral density defined as 

/+oo 
e iut (F{t)F) E dt > 0, (5) 
-oo 

where F{t) is an environmental observable evolving under the Heisenberg evolution governed by the environmental 
Hamiltonian. In the next step of construction we decompose the time-dependent (in the Heisenberg picture) coupling 
operator 5* into the Fourier components 

S(t) = U*(t)SU(t) = Y.Y1 S q {u))e-^ 0+q ^ t 1 n = — (6) 

q£Z {Q} T 

where U(t) is a unitary propagator defined by the time-ordered exponential 

U(t) = Tcxpj-i^ H(s)d s y (7) 

The detailed structure of is a consequence of the Floquet theory. Namely, one can uniquely define the averaged 
Hamiltonian and its spectral decomposition, 

U(r) = e- iS T, H = J2tk\k)(k\ (8) 
fc 

The differences of quasienergies {ik — £;} form a set of Bohr quasifrequencies denoted by {to} and the operators S q (ui) 
satisfy the relation 

[H, S q (Cj)} = ~CjS q (6j) , S- q (-u>) = S\{Q). (9) 

Those transition operators called also Lindblad operators describe the transitions between the levels of the averaged 
Hamiltonian due to the interaction with environment and assisted by the energy exchange q£l with an external source 
of periodic perturbations. They are the main ingredients of the following LGKS generator describing the irreversible 
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evolution of an open system in the interaction picture. This generator is a sum of local generators each of the LGKS 
type 

^ = E E ^ (io) 

q£Z {ui>0} 

C qaP = i{G(w + gfi)([S,(w)p I Sj(^ + [St(p),pS q (Q)])} (11) 

Notice that this construction can be easily generalized to more complicated interaction Hamiltonians of the form 
Hint = T]j S J ® Fj. Moreover, in the case when an environment consists of independent subsystems labeled by j, the 
construction is additive, leading to the generator 

C = E E E C U with G H H- CP(u) , S q {Co) ^ S((Q). (12) 

3 96Z {O>0} 

Finally, we can write down the Markovian master equation in the Schrodinger picture 

Ml = -i[H(t),p(t)} + C(t)p(t) (13) 

where 

C(t)p = W(i)£(W _1 (t)/3) , U(t)p = U{t)pU\t) (14) 

The properties {?]), (j9)), (| IQf) . and (|T2|) imply that the solution of the master equation fp~3|) can be written in a factorized 
form, 

p(t)=W(i)e t£ p(0). (15) 

The decomposition (|15|) is very useful, because it allows us to find the long-time evolution of the system. Namely, if 
C possesses a unique stationary state p (i.e., Cp = 0), then for t —> oo, e tc p(0) — > p and the system asymptotically 
reaches the periodic limit cycle (which for some models is degenerated to a single stationary state) described by the 
formula 



p(t)=U{t)p = p(t + T). (16) 

Periodicity follows from the commutation of p with U (t). In the steady-state regime, only the properties of interaction 
picture generators C? qa matter. In particular, those generators preserve the diagonal (in H basis) elements of the 
density matrix, which is a consequence of (|9]). Therefore, the calculations of basic quantities for the steady state 
involve only the corresponding "classical" Pauli master equations. This property simplifies the analysis of various 
models and often allows for analytical solutions. 



III. THERMODYNAMIC AL PROPERTIES 



The rich structure and nontrivial mathematical properties of the master equations discussed above allow us to 
define consistently heat currents and stationary power, and derive rigorously the second law of thermodynamics for 
driven open quantum systems. 



A. Heat currents and power 



The crucial property of the master equation (fTS"]) is its additivity with respect to the independent system-reservoir 
couplings labeled by j, the Bohr quasifrequencies {ui > 0}, and the harmonics {q£l} of the periodic perturbation. 
The similar additivity is expected for the heat current flowing through the system, which should be a sum of local 
time- dependent heat currents J qSl ify. They are defined as 

Jm = ^-^TT[(C qi (t)p(t))H(t)] (17) 
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where 

H{t) = U\t)HU{t). (18) 

The formula (|1T[) . albeit a bit complicated, possesses a clear physical interpretation. The temporal rate of the energy 
change entirely due to the j-th component of the bath, g-th harmonic of the external periodic perturbation, and the 
transitions between the levels of the temporal averaged Hamiltonian H(t) separated by the Bohr quasifrcquency Hi is 
given by Tr \iC qQ {t)p{t))H{t)\ . The "renormalizing" coefficient (Lj + qfl)/ui takes into account that the system energy 
change u is accompanied by the energy exchange qQ with an external source of the periodic perturbation. 

We assume now that the j-th reservoir is a thermal bath at the equilibrium state characterized by the temperature 
T- 7 and the Kubo-Martin-Schwinger condition for an arbitrary choice of the bath observable F in ([5]) 

G'\-uj) = e~" /T3 G(uj). (19) 

The relation (fT!))) implies the existence of Gibbs-like time-dependent stationary state p J qG j(t) for any local time- 
dependent Schrodingcr picture generator C J qa (t) and similarly for a constant interaction picture one, 

4a(*)#o(*) = °> 4^ = ° ( 2 °) 

exp{-^±$M} 
Trexp{-^Wj 

The the local current then can be rewritten in the following form: 

Jh(f) = -T j Tr[(C q0 (t)p(t)) ln^(i)] . (22) 
and the heat current flowing from the j-th bath is defined as 

w) = £ E 4>(*)- ( 23 ) 

qez {Q>0} 



B. The second law 

The most important result which follows from the theory presented above is the second law of thermodynamics 
written as a positive entropy production, 

|s(p(*))-£^'(*)>0 (24) 

3 

where S(p) = — Tr(plnp) and the inequality ([24]) is a consequence of f|22f) and the Spohn inequality 
Tr[(£p)(ln p — In p] ) < is valid for any LGKS generator C and its stationary state p [![. 

We are often interested in the steady-state (limit cycle) properties. In this case, the formulas simplify, whereby 
the local heat currents and the entropy become time independent. Using the following notation for a stationary local 
current: 

p ^± TT[(cum (25) 

— — UJ 

q<£Z {ui>0} 

one can write the second law for the stationary regime as 

Efj^°- (26) 

3 
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Due to the energy conservation, one can define the stationary power of the external periodic force executed on the 
system as 



v = -^j j . ( 27 ) 



j 

One should note that due to the properties © , (jTTj) , and ([2~Tj) , the interaction picture local generators C 3 qa transform 
the diagonal part [in H basis ([5)1] of the density matrix into a diagonal one and all stationary states p 3 q ^,p are also 
diagonal, which drastically simplifies all computations for particular examples. 

IV. EXAMPLES 

We consider models of a TLS system with the Hamiltonian © coupled to (i) a single zero-temperature bath 
(electromagnetic field) and (ii) two heat baths at different temperatures. The first model provides a description of 
fluorescence in a vacuum and will be used to compare our approach with the standard theory based on the Eqs. 
(P)-®. The second one is expected to describe a microscopic heat pump powered by laser radiation. 

A. Preliminaries: Transition operators for couplings between TLS and environment 

The propagator ([7]) for the Hamiltonian ([2]) can be written as a product of two one-parameter unitary groups, 

U(t) = Ui{t)U 2 {t), Ui{t) = cxp (-^itQa 3 ^ , U 2 (t) = exp { -it Qact 3 + ga 1 ^ j . (28) 

In order to prove ([2"5j) it is enough to differentiate two matrix- valued functions-£/(i) given by ([7]) and Ui(t)U 2 (t)-to 
show that they satisfy the same differential equation. The mean Hamiltonian H reads 

H=±Aa 3 +ga 1 , (29) 

where A = ujq — f2 is the detuning parameter. Next we diagonalizc this Hamiltonian, such that H = ei \4>i) ((f>i \ + 
£2 102) (02 1, with 

ei = -I(0-n fl ), e 3 = ~(Q + n R ), (30a) 
1 ( A + Q R \ . , . _ 1 / A- fl R 



^' Twrtm y * )• * '■ (30b) 



where flu = y 7 4g 2 + A 2 is the Rabi frequency (we assume always SIr < CI). Obviously, the vectors {<pi, 02} form an 
orthonormal basis and, from now on, all matrix representations of appearing operators shall be explicitly written in 
terms of this basis. We choose the following convention for those basis vectors: 

lfc>=(j)> !&>=(?)• (31) 
First, we give the detailed structure of transition operators for two types of TLS-environment couplings. 



1. Coupling by a 1 ® F interaction Hamiltonian 
We compute now explicitly the Fourier decomposition ^ of cr 1 (i). Using ([2"5)l we obtain 
a 1 {t)= ( -7T— cos fit cos Qnt — sin fit sin £Iri\ a 1 — f cos fl^t sin fit + — — cos fit sin flat j a 2 + ( cos fit J er 3 , (32) 

\flR J \ "fl / \Ur / 
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and then its matrix representation reads [in the new basis (1311) ] 

cos ilt + iil r sin Sit) 
° [ >~ fT^ V e ~ mRt ( A cos nt ~ in R sin nt ) - 2 .9 cos m 

By inserting exponential functions into trigonometric ones, one can find all components of the Fourier decomposition 
of this operator. The only positive frequencies correspond to Q = 0, ±S1_r and q — 1 and yield the transition operators 
of the form 

with the symmetry S' 1 (— w) = S' 1 (a;)T [compare cq.©]. 




2. Coupling by o A ® F interaction Hamiltonian 

Just as previously, we have to perform all the necessary mathematical steps once again, this time for the matrix er 3 
to obtain the corresponding transition operators. Again, 

o z {t) = (-^cosClRtj a 1 + (^smQ R t\ a 2 + ^-a 3 , (35) 



with the matrix representation 



3/^ 1 ( A -2ge iQRt , 



Notice that now only two non-negative frequencies appear, and the Rabi frequency £Ir, and with corresponding 
transition operators, 



, 3 , n . _ A / 1 o3fn , 2 9 ( 



and S 3 (-Sl R ) = S' 3 (^i?) t 



B. Example 1: Resonance Fluorescence in vacuum with detuning 

As a first example, we will consider the case of fluorescence in vacuum [lj. A TLS with periodic driving of 
frequency f2 described by the time-dependent Hamiltonian @ is weakly coupled to an electromagnetic bath at zero 
temperature (vacuum). The interaction Hamiltonian between the TLS and the electromagnetic field is given under a 
dipole approximation as 

H int = a 1 ®(ot(/) + o(/)). (38) 
Notice that the usual rotating wave approximation (RWA) was not performed. The short-hand notation 

a ^)= E /OTMk)d 3 fc, (39) 



where k is a wave vector, \i denotes polarization, and /^(k) is a properly chosen form factor used with the electro- 
magnetic field operators satisfying [a M (k), o M 'T(k')] = <5 MA1 '<5(k — k'). The Hamiltonian of electromagnetic bath is given 
by 



h r =Y j /-(kKt (kK(k)d 3 fci w (k) = c|k|. 

M=±1 R3 



(40) 
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The standard choice of the form factor / M (k) = z(27r) 2 ^/c|k|/eo(d • e^), where d is the transition dipole and et M 
determines polarization, yields the spectral density at vacuum state 

G(u>)=Au>\ A=-^. (41) 
For convenience, let us introduce some useful notation, 



5 ° = {it) Gm > 6± = {^r) Gin±n ^ (42) 

71 =S-+6 + , 72 = i(«J_ + 5+ +6 ). (43) 

Performing all the necessary calculations and summing up all parts like in the prescription in (|10[) . (|lip. and (|12l) . 
one obtains the resulting semigroup generator £ acting on the reduced density matrix p 1 (t) of TLS in the interaction 
picture and, finally, the associated master equation in matrix elements [pii(i) + P22W = 1] 

= .g+fat) + ^ (t)> = _ 72 ^ (t) (44) 

with the solution 

p( 1 (0= e - 7lt pii(0) + — (1-e-^ 4 ), p{ 2 (t) = e~^ Pl2 (0) (45) 

7l 

Using those results, we can calculate the fluorescence power spectrum. We apply the standard formula describing 
the power spectral density function of incident light scattered by a TLS coupled to an electromagnetic reservoir, which 
can be found, for example, in Ref. [l[ 



00 

J(w) ocRc J e- luJt tY{a + U(t)e tc {<j-p)} dt, (46) 


where a ± = \(o l ±icr 2 ), U{t)e tc is a quantum dynamical map in the Schrodinger picture and p stands for stationary 
state of semigroup generator, satisfying Cp = 0. After straightforward calculations, one can find the formula 

Kfj) oc ~** + ) x <K W - n) (47) 

, , / Ag 2 5+ \ ( 4 5 2 <5_ \ 1 7l 



. (n« - A) 2 7l y V (n B + A) 2 7l 7 7T 71 2 + (w - O) 2 
45 2 <S+ 1 72 



x 



(0^ + A) 2 7l 7T 72 2 + [« - (fi - Ok)] 



2 



,2 



5_ 1 _ 72 



(Ofl-A) 2 7l 7r 7 2 2 + [ W -(0 + O i? )] 2 ' 

which clearly possesses the structure of the usual Mollow spectrum {Mollow triplet) , with Lorentzian peaks centered 
around frequencies fi (central peak) and fi ± FIr (side peaks) . The Dirac delta contribution originates from possible 
clastic scattering of incident light without changing the energy of scattered photons. All terms are products of 
Lorentzians (or Dirac delta) with normalized integrals and weights describing the relative contributions to the total 
power of the scattered radiation. Although, qualitatively, the obtained line-shape is standard, the details of the 
coefficients differ from those found in the literature, in particular computed using (|3|) (see, e.g., Ref. [l|). Our 
formulas contain more details of the model including, for example, the frequency dependence of the decay rates. The 
plot of the obtained power spectral density function (|4"T)) is presented in Fig. [T] Notice that the chosen ratio g/fl and 
the value of A are very high in comparison to known physical examples in order to make visible the asymmetry of 
the line shape. 
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FIG. 1: Power spectral density function (|47[) in a process of nonresonant fluorescence in vacuum (fi = 0.85, A = O.Of , g — 0.075). 
The chosen ratio g/Q. and the value of A are very high in comparison to known physical examples in order to make visible the 
asymmetry of the line shape. 



C. Example 2: A microscopic heat pump driven by a laser 

As a second example, we will consider a model of a TLS, periodically driven by a laser field and weakly coupled to 
two heat baths: an electromagnetic (photon) bath at the temperature T e > and a dephasing bath at the temperature 
Td > 0. The TLS-bath interaction Hamiltonians differ for both cases: The photon bath is coupled by the Hamiltonian 
of the form (|38|) (i.e., er 1 g) F -coupling) while the coupling to the dephasing bath is of the er 3 g) F -type. The second 
form of coupling means that without external driving the second bath cannot cause transitions between the TLS 
energy levels and produces only dephasing (pure decoherence) of the TLS state. This is typically due to the energy 
mismatch between high energy level splitting of the TLS and low energy excitations of the bath. The external periodic 
driving together with the presence of the electromagnetic bath generate energy transfer between the laser field and 
both heat baths mediated by the TLS. It follows from the formula (|3T>)) which shows that the driven TLS is coupled 
to the dephasing bath through oscillations with Rabi frequency which is much lower and, hence, can match bath 
excitations. 

The following examples of physical systems satisfy the assumptions of this model: 

(1) An optically active atom immersed in a dense buffer gas. The collisions with the buffer gas are elastic; hence, 
the condition of a pure dephasing bath is satisfied. Recently, both cooling (for red-detuned laser radiation) and 
heating (for blue-detuned laser radiation) have been observed using rubidium atoms in argon buffer gas |{|. 

(2) A quantum dot interacting with acoustic phonons. As the optical frequency of the dot is much higher than the 
Dcbye frequency, a "bottleneck effect" appears [l(|, which suppresses the energy exchange between the dot and the 
phonon bath. Again, a proper detuning of the laser beam should generate cooling or heating of the phonon bath. 

(3) A metallic nanoparticle supporting a TLS and coupled to a macromolecule. The molecule vibrational degrees 
of freedom form a dephasing baths. Recently, such hybrid systems driven by laser radiation were investigated in the 
context of molecular electronics [TTj . It seems that such systems could also be used as nanomachines (heat pumps or 
engines). 

The quantitative analysis of the discussed model is based on the master equation and the thermodynamical for- 
malism developed in Sees. II and III. Similarly to the previous example, we can use the formulas for transition maps 
and construct the interaction picture generator. We use also a short-hand notation similar to that from the previous 
section, 

5o = (^) 2 G d (n R ), s± = (^|±^) 2 G e (tt±), n ± = n±n R (48) 

where G e (oj) and G d {oj) are spectral densities for the electromagnetic bath and dephasing bath, respectively. 

To obtain the expressions for steady-state heat currents we need only the following equation for the diagonal matrix 
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elements in the interaction picture (p 22 — 1 — p\\ ) 

d P [S) 



= -[So + e —5-+S + )p I 11 {t)+(e ^S + 8^ + e — 8+)p I 22 (t), 



(49) 



dt 

Finally, we compute heat currents for both baths using the formulas (|25|) and the stationary state obtained from ((49 
The matrix elements of p satisfy 



P22 



e T d § + S- + e T = 8 + 



S + e T eL| 5 + 

The formulas for the heat currents flowing from the dephasing and the electromagnetic baths read 

K 



1-e 



8+ e r «™ - e T c, 



J' 



1 



QrK, 



2 fe~ — - 1J 



i mu, + e 



1 5 ^- 



e ^ ] n + 6 8 + 



where 



> 0. 



(50) 

(51) 
(52) 

(53) 



The expressions (|51j) -(|53 |) can be used to analyze particular experimental settings. Here, we show only analytically 
the phenomenon of switching from the cooling to the heating regime observed in the experiment Q . We begin with 
the simplifying assumptions which are typically valid in quantum optics domain, 



n > T e , r d , n R « n, n R < T d . 



(54) 



The last inequality holds for small detuning A = luq — only, but we are interested in the region where detuning 
changes its sign. Under the conditions (|54| we can put 



and using ([5"l"1) - (l5"3]) we obtain 



while 



0, 



V A, V = 



i, G e (n ± ) ~G e («), 

2S Q + 5- +5+ ' 



J" 



2S + S-+ 8+ 

As a consequence, we obtain two regimes for small detuning |A|: 



(55) 



(56) 



(57) 



(a) Cooling regime: For red detuning (i.e., 51 < u>o), A > and J d > 0, which means that heat flows from the 
dephasing bath to TLS. The system acts as a heat pump cooling the dephasing bath and heating the electromagnetic 
one at the expense of the work provided by the laser field. 

(b) Heating regime: For blue detuning (i.e., f2 > luq), A < 0, both heat currents J d ,J e < 0, which means that 
both baths are heated at the expense of the laser field. 



V. CONCLUSIONS 



The presented formalism based on Markovian master equations describing periodically driven quantum open systems 
weakly coupled to heat baths is very flexible and possesses a large number of possible applications. It is consistent 
with thermodynamics and leads to simple forms of kinetic equations, heat currents, and stationary power. This 
allows us to study, to a large extent analytically, various designs of microscopic machines, such as refrigerators and 
engines. The model briefly discussed here as Example 2 will be used in forthcoming papers for a detailed quantitative 
explanation of the experimental results of @ and as a basic ingredient of a solar engine model. 
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